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Abstract. We compare the ultra-violet, optical, and far- 
infrared emission for a sample of 135 spiral galaxies in order 
to address the widely debated problem concerning the opac- 
ity of spiral disks. We find that the re-radiation of the dust, 
estimated from the far-infrared emission, is on average only 
31±1 percent of the bolometric luminosity of a spiral galaxy, 
indicating that less than one third of the stellar radiation is 
absorbed and then re-radiated by dust in a spiral disk. Ap- 
plying a radiation transfer model which assumes a 'Sandwich' 
configuration for the spiral disk, and fully takes into account 
the effect of scattering, we find for our sample a median of the 
face-on blue (4400A) optical depth tb = 0.49 and the mean 
< Tb >= 0.60 ± 0.04, indicating that most spiral galaxies in 
our sample are not opaque for blue light (tb < 1). 
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1. Introduction 

The question whether the disks of spiral galaxies, which con- 
sist of stars, gas and dust, are transparent (optically thin) or 
opaque (optically thick) in optical wavelength range, is cur- 
rently a hotly debated one. It is equivalent to the question "Are 
the spirals intrinsically very different from what we see from 
the Earth?" Our understanding of the nature of these astro- 
nomical systems, which contain most of the visible material of 
the universe, depends fundamentally on the ultimate answer to 
this question. The classical approach to the problem is based on 
the statistical studies of the dependence of the optical surface 
brightness /j, on the view angle i (i = for seen face on), under 
the hypothesis that for a transparent disk /j, depends sensitively 
on i, while for an opaque disk it does not (Holmberg 1958). Ap- 
plying this method to a sample of more than 9000 disk galaxies, 
the largest sample to date, Valentijn (1990) claims that most 
spiral disks are opaque (tb > 1), and therefore are intrinsically 
much brighter than observed. However, this approach is rather 
indirect and is influenced sensitively by biases in the sample 
selection (Disney 1992). Consequently, conflicting results have 
been obtained by different authors (Holmberg 1958; de Vau- 
couleurs 1959; TuUy 1972; Davis et al. 1989; Valentijn 1990; 
Burnstein et al. 1991; Boselli &; Gavazzi 1994). 

Modern astronomy has opened the entire domain of the 
electromagnetic waves, from radio (A > 1 cm) to 7-ray 
(A < O.IA). All the energy radiated by a spiral disk can, in 
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principle, be recorded. The problem of the optical thickness 
can be addressed directly by comparing the escaped emission 
of stars, the major sources of emission of galactic disks, with 
the re-radiated emission of dust which causes the extinction. If 
the disks are indeed opaque, the latter should be much larger 
than the former, and vice versa. Stars emit predominantly in 
ultra-violet (UV: 500-3650A), optical (3650-9000A), and near- 
infrared (NIR: 9000 A - 3//m); dust emits in mid- infrared (MIR: 
5-40//m) and far-infrared (FIR: 40-1000//m). In this paper, we 
compare these two emissions for a sample of spiral galaxies 
and find that most spiral disks in our sample are not opaque 
{tb < 1). 

2. Dust emission to bolometric emission ratio 

Our sample is selected from nearby spiral galaxies detected in 
two recent UV experiments: FAUST (Bowyer et al. 1993; De- 
harveng et al. 1994) is a far-UV telescope (central wavelength 
1650A) flown in 1992 on board of the space shuttle Atlantis; 
and SCAP (Donas et al. 1987) is a balloon-boarded UV ex- 
periment (central wavelength 2030A) carried out in Marseille 
during the early 1980s. Galaxies are included if passing the 
following two criteria: (1) detected by IRAS in the 60 and 100 
/im bands, (2) morphological types from Sa to Sm. The sample 
contains 135 galaxies. The monochromatic UV fluxes are taken 
at 2030 A, the FAUST data are translated to this wavelength 
according to the calibrations of Deharveng et al (1994) (a sys- 
tematic difference of 29% is found between the FAUST and 
SCAP fluxes). The IRAS 60//m and 100//m fluxes are taken 
from the IRAS data base. The monochromatic 4400A fluxes 
f440o are taken from the Third Reference Catalog of Bright 
Galaxies (de Vaucouleurs et al. 1990, hereafter RC3). Both the 
UV and the blue fluxes are corrected for the foreground Galac- 
tic extinction using the blue extinction tabulated in the RC3 
and the extinction law of Savage and Mathis (1979). 

In order to estimate the integrated UV flux fuv (912 - 
3650A) from the monochromatic flux f2030, the average broad 
band UV spectra (in the range of 1400 - 3650 A) reported 
by Pence (1976) for Sa-Sab galaxies, and by Coleman et al. 
(1980) for Sb-Sbc, Sc-Scd, and Sd-Sm galaxies are extrapo- 
lated to the range of 912 - 3650A, and then are adopted for the 
corresponding types. The optical-NlR fluxes fop-nir ^'■^ S^H" 
erally estimated from the blue fluxes using the average broad 
band optical-NIR spectra for different Hubble types estimated 
from the unweighted mean of the empirical spectra collected 
by Yoshii &; Takahara (1988), except for 85 galaxies whose 
B — V colors are taken from RC3, and for 52 galaxies whose 
U — B colors are taken from the same source. In a first time we 



neglect the ionizing UV radiation (< 912 A) which is predomi- 
nantly absorbed by dust. Finally, we estimate the total escaped 
radiation of stars from a disk galaxy by fstar = fuv + fop-nir' 
The dust emission is estimated from IRAS fluxes feo/i and 
fioo/i, namely the radiation in the wavelength range of 40 — 
120//m. An important issue is the radiation of the cold dust 
grains (~ 15K) at wavelengths > 120//m. There is evidence 
from IRAS studies and mm/sub-mm observations that such 
cold grains dominate the total mass of dust (Devereux &; Young 
1990; Kwan fc Xie 1992; Bothun fc Rogers 1992; Chini fc 
Kriigel 1993). However, because of the strong dependence of 
FIR luminosity on grain temperature (L oc T , n=l ~ 2), 
the contribution of cold grains to the total emission of the dust 
is not very significant (Devereux &; Young 1990; Devereux &; 
Young 1994). For example, recently Guelin et al. (1993) find 
that the 1.3mm continuum flux of NGC 891 is about a fac- 
tor of 8 more than a simple extrapolation of the IRAS fluxes 
suggesting the existence of significant amount of cold dust, 
but the luminosity of this cold dust (Led) is only 1/3 of that 
of the 'warm dust' (Lwd). Nonetheless, large uncertainties are 
present when only IRAS data and 1mm data are compared 
(Devereux &; Young 1990; Valentijn 1991), and different con- 
clusions have been reached by different authors; for instance 
Chini et al. (1986), fitting the IRAS data and their 1.3mm 
data using a two-component model, found that the cold dust 
emission (18K) dominates the warm dust (53K) emission. The 
best way to estimate the cold dust radiation beyond 120//m 
is perhaps to compare the IRAS data and the sub-mm data 
(mostly in 300-800//m) which are much closer to the peak of 
the dust emission (likely in 100-200//m). Recently, using avail- 
able sub-mm data (Chini et al. 1986; Stark et al. 1989; Eales 
et al. 1989), Kwan fc Xie (1992) estimated the total FIR lumi- 
nosity (40-1000//m) for 11 late-type galaxies. Adding NGC 891 
(Guelin et al. 1993) and the Galaxy, the FIR-sub-mm spectrum 
compiled by Mazzei et al. (1992, see also Chini &; Kriigel 1993), 
we find a mean of L(40-1000//m)/L(40-120//m) = 1.39 ± 0.21 
for these 13 galaxies. 



Fig. 1. Diagram of fdust/fflr ratio versus Iog(f6o,i/fio0/i) for: 1) 11 
galaxies with sub-mm data taken from Kwan & Xie (crosses); 2) 
NGC 891 (open square); 3) the Galaxy (solid square). 

There is also 'hot dust' consisting of small-grains/large- 
molecules heated transiently by single UV or optical photons 
to T > a few hundred K. It emits at wavelengths shorter than 
40//m (MIR), and may account for 20-30% of the total dust 
emission (Puget &; Leger 1989; Desert et al. 1990). Most of this 
emission is included in the IRAS 12//m band (8-15//m) and 



the 25//m band (18-30//m). Hence, again for the 13 galaxies 
mentioned above, we estimate the MIR fluxes (8-40//m) using 
the formula 
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where fi2,i and f25^ are in Jy, c is the speed of light, Ai = 12//m 
and 5X1 = T/im, and A2 = 25//m and 8X2 = 25//m. The value 
of 5X2 is larger than the real bandpass of IRAS 25//m channel 
(11.15//m) in order to cover the wavelength range 15-40 /im. 
The total fdust/fflr ratio, fflr being the integrated FIR fluxes in 
the wavelength range 40-120//m estimated from feo/i and fioo/i 
(Helou et al. 1988), is estimated as 



fdusl ^ fmir f(40-1000//m) 

fflr fflr fflr 



(2) 



In Fig.l we plot the fdust/fflr ratio versus log(f6o,i/fioo/i) for 
the 13 galaxies with sub-mm data. The ratio spans from 1.4 
to 2.2, with a mean 1.79±0.25. A strong linear anti-correlation 
between fdust/fflr and log(f60/i/fioo/i) is found (linear correla- 
tion coefficient -0.76). This is easy to understand because the 
warmer the FIR color ratio (feo/i/fioo/i), the more contribution 
from the warm component to the total dust emission (Xu &; 
De Zotti 1989; Kwan &; Xie 1992). The dotted line is the linear 
regression 

fdust/fflr = 1.45 (±0.08) - 1.06 (±0.24) x log(f60;./fioo,.) (3) 

We shall adopt this empirical relation for calculations of fdust 
from fflr and feo/i/fioo/i. It should be noticed that this relation 
is derived from a small sample of galaxies, and there might be 
some systematic uncertainties with their sub-mm data (Chini 
&; Kriigel 1992). Our results might be affected by these uncer- 
tainties. 



Fig. 2. Histogram of dust— to— bolometric emission ratio. The bolo- 
metric flux fboi = fstar + fdust- 



Figure 2 is a histogram of the fdust /fboi distribution for 
our sample, where the bolometric flux fboi = fstar ± fdust. The 
sample mean is 0.31±0.01, indicating that on average the dust 
absorbs less than one third of the stellar radiation. This is 
consistent with the result of Cox and Mezger (1989), who find 
that the dust emission is only one third to one fourth of the 
bolometric luminosity of the Galaxy. 



3. A detailed radiation transfer model 

In the literature, the opacity of disks is usually studied in terms 
of tb , the optical depth for blue (4400A) radiation. In the fol- 
lowing we estimate tb for the galaxies in our sample using a 
detailed model for radiation transfer and dust heating. 

We consider three contributions to the heating of dust: 1) 
dust heating due to the ionizing radiation (A < 912A), 2) due to 
the non-ionizing UV radiation (912A < A < 3650A) and 3) due 
to the optical radiation (3650A < A < 9000A). We assume that 
the ionizing radiation mainly heats the dust around the massive 
stars with which the Hll regions are associated. Subtracting 
this contribution from the total dust radiation, we estimate the 
optical depth of the disk by comparing the remaining (diffuse) 
dust emission, the non-ionizing UV and the optical radiation 
for a sample of galaxies using a radiation transfer model. 



Fig. 3. Schematic illustration of definitions of several parameters 
used in the radiation transfer model, e denotes the volume emissivity, 
K the volume absorptivity, ^' = cosi', and ri = / Kdz 

I/O 

The ionizing UV radiation (A < 912A) can be estimated 
from the Ha emission. For 34 galaxies in the sample, the Ha 
fluxes fnct are available from Kennicutt and Kent (1983) and 
Romanishin (1990) (by order of preference). After correcting 
isa for 0.8 mag internal extinction (Kennicutt 1989), the ion- 
izing flux is estimated as follows: 



fLyc = N(Lyc)x < E(Lyc) >= 33.9 x fnc 



(4) 



where N(Lyc) = 7.19 10^V0-'''5 x fac is the number flux (in 
units of cm~ see" ) of Lyman-continuum photons (Mezger 
1978; Lequeux 1980), and < E(Lyc) >= 2.13 x E(Lya) = 
3.49 10~^^ erg is the average energy of Lyman-continuum pho- 
tons (Mezger 1978). Assuming that 80 % of Lyman-continuum 
emission is absorbed by dust, directly or indirectly (via emis- 
sion lines) and in or outside Hll regions, the contribution from 
the ionizing UV radiation to fdust can be estimated: 



fiyc 



0.8 X fLyc 



(5) 



For 11 early spirals (Sa-Sbc) with fna, we find that on average 
the ionizing UV contributes 14±2% of total fflr. For the 23 later 
galaxies the contribution is 20±1%. We estimate the values of 
^du'st ^°^ galaxies without Ha data using these means. The rest 
of fdust , fi = fdust — fdu'st ' ^^ ^^^ ^° non-ionizing UV and optical 
radiation. 

The radiation transfer model for the non-ionizing UV and 
optical radiation of disk galaxies is illustrated in Fig. 3. It is 
based upon the model of Xu &; Helou (1994) which has been 
applied to the heating of the diffuse interstellar dust in M31. 



The model takes the effect of scattering fully into account in the 
sense that scattered light of any order has been calculated using 
an iteration procedure from lower order scattered light (van de 
Hulst &; de Jong 1969). An infinite-plane-parallel geometry is 
adopted for the radiation transfer problem, and the radiation 
sources (stars) and dust are assumed to be smoothly distrib- 
uted. We allow for different thicknesses of the star layer and of 
the dust layer ('Sandwich model'). The model calculates for a 
given wavelength, at a given point in the disk and in a given 
direction, the radiation intensity l;\(ri,// ) as a function of the 
radiation source function So, a and the optical thickness of the 
disk T\. From this we predict the ratio between the radiation 
absorbed by dust (dust-heating) and the radiation escaping the 
disk and eventually being observed, as a function of the optical 
depth of the disk and the view angle: 



Qx{tx,Ii) 
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where 1;\(0,//) + Ei/// is the prediction for the light observed 
(see Fig. 3 for the definition of Ei), /j, = cosi the inclination, t\ 
the face-on optical depth, w the solid angle, and a the albedo. 
The factor (1 — a) gives the ratio between the absorption cross 
section and the extinction cross section. 

We split the dust heating into that due to the non-ionizing 
UV radiation (UV-heating) and that due to the optical radia- 
tion (optical-heating). The UV-heating is calculated as 



=f2030 X Q(r2030,/i)x 
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where fdust i^ the dust radiation intensity due to the UV heat- 
ing; f2030 the 2030A flux, and r203o the optical depth at 2030A; 
(m(A) — m(2030A)) is estimated from the adopted UV spec- 
trum. Similarly, the optical heating is 
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The optical depth at different wavelengths (including r203o ) are 
related to tb through the adopted extinction curve. 
The settings of the dust-heating program are 

1) Sources of non-ionizing UV radiation (i.e. B stars) have the 
same scale height as dust (Ei = E2 =0), while sources of 
optical radiation (i.e. old stars) are in a disk twice as thick 
as the dust layer so that half of these sources lie outside 
the dust layer, i.e. Ei = E2 = 0.5So,a x t\. 

2) The UV albedo = 0.18 (Hurwitz et al. 1991), and the scat- 
tering asymmetry factor g = (isotropic scattering). 

3) The optical albedo and phase function are taken from 
Mathis et al. (1983). 

4) The UV extinction curves are taken from Pei (1992) for 
spiral galaxies. 

5) The optical extinction curve is taken from Mathis et al. 
(1983). 

For each galaxy in our sample, the face-on optical depth tb 
is estimated from the equation 

fdust = f^st + [C^t {tb ) + ill, irs)] (9) 

where i/^^, is the dust emission due to ionizing UV flux. The 
results are plotted in Fig. 4. The histogram of the tb distribu- 
tion for all galaxies in the sample has a median value of 0.49, 
indicating that most of them are not opaque (tb < 1). The 
mean value < tb >= 0.60 ±0.05 is 'robust', i.e. the 3(7+ points 
are not included in the calculation. The distribution is rather 



Fig. 4. Histograms of tq (face-on optical depth for radiation at 
4400A) distribution. Hatched: late spirals (Sc— Sm). 



asymmetric with a sharp peak at tb ~ 0.2 and a long tail at the 
large tb side. Therefore the mean is larger than the median. 
The histogram (hatched) of the tb distribution for late-type 
(Sc-Sm) spirals is similar to that of the total sample, with the 
median= 0.36 and the mean= 0.43 ± 0.04. Early type spirals 
(Sa-Sbc) are less transparent, their tb distribution is flatter 
than that of the late types with a median value of 0.86 and a 
mean equal to 0.79 ± 0.08. It is interesting to note that for 59 
Sb-Sc galaxies in our sample, we find a median tb = 0.90 and 
< Tb >= 0.95 ± 0.07, slightly less than the result of Valentijn 
(1991) for the same types of galaxies (tb ~ 1.3). 

Since our sample is basically UV selected, there might be 
a bias in the sense that less extinguished galaxies are more 
likely to be included. However, we find for our sample a mean 
Lflr/LB = 0.79 ± 0.10, where Lflr is the integrated FIR lumi- 
nosity in the wavelength range of 40-120//m, and Lb the blue 
luminosity (i'L„(4400A)). This is consistent with the FIR-to- 
blue luminosity ratio of UGC galaxies (< Lflr/LB >= 0.72, 
Bothun et al. 1989), indicating that the bias, if existing, is not 
serious. A more detailed discussion of the extinction at UV 
wavelengths derived from this model is postponed to a subse- 
quent paper (Buat &; Xu, in preparation). 

4. Discussion 

Disney et al. (1989) exploit a similar method to the problem of 
the optical thickness of the spiral disks as we do in this letter. 
However, they compare the FIR emission of spiral galaxies only 
with the optical emission, neglecting the energy absorbed by 
dust from the UV radiation. They concluded that the data are 
consistent with the spiral disks being generally optically thick 
(tb > 1). We reach a different conclusion because of includ- 
ing the contribution of the UV radiation to the dust heating. 
Indeed, we find that the UV radiation dominates the optical 
emission in dust heating: the fraction of dust heating due to 
non-ionizing UV (92lA < A < 3650A) is 60±9%, and that due 
to optical radiation (3650A < A < 9000A) is only 21 ± 8%. 
The rest is due to the ionizing radiation (A < 912A). This is 
consistent with the result of Xu (1990) who found that the 
non-ionizing UV radiation can account for 57-76% of the ra- 
diation of the cool (diffuse) dust in spiral disks. The heating 
of dust by the non-ionizing UV radiation, radiated by interme- 
diate massive stars (~ 5M0 ) plays an important role in the 
remarkably tight and universal correlation between the FIR 
and radio continuum emission of spiral galaxies, because the 



same stellar population may be responsible, through supernova 
explosion, for the acceleration of the cosmic-ray electrons ra- 
diating the nonthermal synchrotron radiation of these galaxies 
(Xu et al. 1994). 
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